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ABSTRACT 

We analyze the primordial protogalactic gas to investigate the role of 
thermal instabilities in the formation of condensations that could contribute 
to the MACHO mass function. We find that the previrialized protogalaxy 
could have undergone a thermally unstable regime which defines a peak in 
mass condensations starting around 0.1 Mq after being reprocessed by the 
environment which may be identified as giving origin to the MACHOs. Effects of 
magnetic fields on this instabihty-driven MACHOs scenario are also addressed. 

Subject headings: Galaxy : Halo : Dark Matter - Stars : Low-Mass, Brown 

Dwarfs 
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1. Introduction 

Two years after the first reports of candidate microlensing events (Alcock et al. 1993, 
Aubourg et al. 1993, Udalski et al. 1993) the field is rapidly growing and bocoming a 
major tool to investigate the morphology and history of our galaxy (see Zylberajch 1994 
and Paczyhski 1995 for reviews). Particularly intriguing is the fact that a massive halo 
composed of MACHOs seems not to be favored by the observations (Alcock et al. 1995, 
Gates et al. 1995). The nature of dark objects so far responsibles for the events towards 
the LMC needs also to be addressed. It is not clear as yet if they belong to a disk of our 
own galaxy, to the galactic halo itself, or even to the LMC (Sahu 1994). We shall adopt 
in this work the hypothesis that MACHOs are indeed located in the halo and attempt to 
outline a scenario that could have leaded to their formation. 

A number of works have stressed the possible relevance of thermal instabilities (Tl) for 
globular cluster formation during the collapse of the protogalactic cloud (e.g.. Fall & Rees 
1985, Vietri & Pesce 1995 and references therein). While these scenarios address an epoch 
in which the gas has already been shock-heated to its virial temperature T ~ 10^ K, we will 
be mainly interested in a previous, optically thin stage of the cloud at T ~ few x 10^ K. 
We also note some previous works dealing with the protogalaxy which have addressed the 
question of Pop III star formation (e.g. Silk 1983). The motivation for a reasessment 
of such scenarios is stronger than ever, particularly now that an initial spectrum can be 
calculated and its reprocessing characterized observationally by the growing body of data. 
To stress this point, let us consider, as in those previous works, a quasi-static, isothermal 
protogalactic cloud with an estimated virial temperature T ~ iJ,v^/2kB, or about 10^ K. It 
is known that at earlier stages the faster cooling at the center would have destroyed the 
dynamical equilibrium and induced the collapse of the cloud, mainly because the equilibrium 
requires the opacity of the gas to be large enough in order to delay radiative losses. On the 
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other hand, an optically thin gas will cause the cloud to cool more quickly, so that at least 
for timescales of the order of the free-fall time tjf ~ {GpY^/'^ the protogalactic gas will 
remain at a temperature much smaller than the virial one. There is also an observational 
motivation to consider a previous epoch of the galaxy; namely that a substantial fraction 
of the halo mass is in the form of field stars, not in globular clusters at a ratio 1/100. 
Furthermore, the most metal-poor clusters have [Fe/H] ~ —2.5, while some field stars 
have been detected with [Fe/H] ~ —4, suggesting that they were formed before globular 
clusters. 



2. Thermal instability 

Motivated by the above discussion we consider a previrialized collapsing cloud 
with mass M = 10"^^ Mq, radius r > 50kpc (corresponding to an average density 
p ~ 10~^^ g cm~^) and T — few x 10^ K. Such a cloud would have started to collapse at 
2; < 10 when its mean density was about 5 times the background density (Rees & Ostriker 
1977). 

At temperatures below T ~ 3 x lO^X the primordial gas of the cloud (which has a 
metal abundance Z ~ 0) coohng through radiative thermal processes (coUisional excitation 
and recombination) is stable to condensation formation by isobaric thermal instability (e.g. 
Field 1965, CorbeUi & Ferrara 1995). We assume that the onset for the thermal instability 
occurs at T ^ 3 x lO^X. We also assume a perfect gas equation of state P — RpT/p , 
where fi =< m > /tjih and p,T are the equilibrium density and temperature of the cloud. 
A perturbed gas element is subject to an isobaric thermal instability growing to the first 
order at a rate (Field 1965, see also Begelman 1990, Gouveia Dal Pino & Opher 1989a,b, 
1990, 1991, 1993) 

n 
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where /i is the mean molecular weight (= 0.59 for an ionized gas of primordial composition), 
7 = 5/3, X = 0.77 X 10-6 is the thermal conductivity term and L = L{p, T) - G{p, T), 
i.e., cooling minus heating, is the net energy loss per unit time and mass, which vanishes 
at equilibrium. Note that, strictly speaking, eq.(l) holds for the short wavelenghts of the 
perturbation. A very important (and uncertain) point is the nature and dependence of 
the heating function G{p,T). The heating function may include the UV flux of an AGN 
phase (likely to have been present in the early galaxy because of the observational evidence 
collected from the bulge) ; turbulence (the Reynolds number Re — 2 rcv/rj is known to 
be ^ 1 for the primordial gas, v being the gas velocity and rj the shear viscosity); and 
MHD waves among the most important sources. Given the uncertainties on the nature 
of the dominant heating mechanisms in the protogalactic gas, we assume as in previous 
works (e.g. Field 1965, Eilek & Caroff 1979, Gouveia Dal Pino & Opher 1989a,1989b, 1990, 
1991, 1993) that the heating function G{p,T) is due to an externally-generated source and 
its dissipation is locally independent of the density and temperature, a condition which is 
known to hold for several of the above mentioned sources. Thus, the derivatives of C in 
eq.(l) will depend only on the cooling function L. 

For a strongly ionized coUisional H + He plasma of primordial composition { X ~ 0.76, 
Y ~ 0.24 and Z ~ 0) at 2.5 x lO^iT <T < lO^K the net cooling is 

L ~ A{T)n]j/p ergg-h-\ (2) 

where uh is the hydrogen number density and the cooling function A(T), evaluated by, e.g.. 
Fall & Rees (1985), can be fitted by 

logA{T) = -22.9 + 2.45 0(T) ea;p(-0(T))2 + ^Ij{T) exp{-^l){T)f + OiT) exp{-e{T)f (3) 

with (t){T) = l.llogT - 4 , il){T) = 2,.7logT - 17.6 and e{T) = lAlogT - 9.2. 
This rather complicated fit can be replaced by the much simpler expression 
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logK{T) = -24.45 + IMlogT - Q.ll{logTY in the restricted region of interest of 
this work 2.5xlO^X < T < IQ^K with sufficient accuracy. Furthermore, since 
\TdC/dT\ <C \pdC/dp\ (see eq.l), we can safely neglect dC/dT in our analysis and adopt a 
constant A(T) ~ 9.6 x 1^'"^^ erg cm^ in eq.(2); this extremely simple form introducing 
errors of at most 10 % in the calculated quantities. 

An isobaric thermal instability grows only if the pressure in the perturbed 
region is equalized in a time-scale small compared to the cooling time, defined as 
tc = {?> kBT){2 jirriH L)"^ . This implies an upper limit for the wavelength of the 
unstable region (e.g. Field 1965) 

A < Xmax - 2TlVst, = 1.1 X 10^"^ {p_25)-' {T^,-rf^cm. (4) 

where Vg — {'^ksT / pmnY^'^ is the speed of sound, p_25 = p / {IQ""^^ g cm~^) , and 
74.7 = 7^/(5 X lO^X). 

On the other hand, it is easy to show that, due to the thermal conduction 
ability to supress the growth of a thermally unstable region, a critical wavenumber 
kmax = '^T^l ^min exists ovcr which n = and it is given by k^ax = [kx {kp — ^t)]^''^ ; where 
kp — ^^rv^t''^'^ ' ~ p. ^-i)K ' ~ ~ — ^Rv^'^^'' ■ Aftsr an elementary manipulation 

with these expressions we find Xmin — '^'^/kmax which is given by 

A^,„ ~ 2.1 X 10^«(p_25)-' (^4.7)'/' cm. (5) 

Since n is small for small k (large A) and vanishes for k > kmax, a k exists for which n 
reaches a maximum. This value is k{nmax) — [kp{kp — /ct)^^:]^^^, or in terms of X{nmax) 

Xin^ax) ^ 1.5 X 10^^(p_25)-' (T4.7)''/' cm. (6) 

The most rapidly growing condensations are those at X{nmax)- They define a typical mass 

Mbiob ~ PbX^irimax) or 

Muob ~ 1.8 X 10-1 (1 + /6e)(p_25)-' in^f'/' Me, (7) 



- 7- 



where fb < 1 is the initial density perturbation relative to the ambient density p. This 
is the most probable mass of the condensations, which is bracketed between the limiting 
values 

MuobiXmin) ~ 4.4 X 10-" (1 + fbe){p-25)~' (74.7)''/' Mo, (8) 

and 

MuobiXma.) - 73 (1 + /,e)(p_25)-' {^.7^^ Mq, (9) 
which define the condensation mass spectrum at early times. 

3. Primordial magnetic field effects 

If a small primordial field is present in the gas B 7^ (which seems unavoidable in 
the pregalactic medium, e.g. Rees 1994) in the direction || to B, the condensation is not 
affected by the presence of the magnetic field and its growth is given by eq. (1). However, 
in the direction _L to B, the thermal conduction K is greatly reduced (because of the e~ 
spiraling between collisions) and the initial collapse of the condensation is more efficient 
(e.g. Field 1965 ; Gouveia Dal Pino & Opher 1989a, 1989b, 1990, 1991, 1993). We can 
therefore determine two scales in this direction X±min and X±max for the instability growth 
with analogous meanings to those of the former section. Defining an effective pressure given 
by the sum of the gas pressure and the magnetic pressure we have (e.g. Field 1965) 

X±rmn ~ 4.1 X 10^^ {B_u)-' {T,,,)"' (1 + ^) ^ cm, (10) 

where va = -6/(4 7rp)^/^ is the Alfven speed, S_i2 = 5/(10""^^ G), and (e.g. Gouveia Dal 
Pino & Opher 1990) 

2 1/2 

A^^„. ~ 1.1 X 10^° (p_25)-^ {T,,,f'' (1 + ^) cm. (11) 

Note that for small B the ratio v\/Vg -C 1 for typical values of protogalactic p 
and T. The masses that bound the spectrum (Ay = X{nmax)) can be estimated as 
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Mbiob{>^±min) ^ Pb'^ >''imin \b ^^^^^ results in 



Muob{X±min) = 4.1 X 10-^ (1 + he){T^.7Y'/^ (B-u)-' (1 + ^) Mq, (12) 



and Mhiob{>^±max) ^ Pb ^^Imax yielding 



MuobiX^ma.) = 31 (1 + he){p-25r' {T^.7f'^ (1 + 4) ^©^ (^3) 



4. Evolution of the blobs 

From eq.(l) we can evaluate the growth time Tg of the condensations in the absence of 

B: 

T, = [l X 10-^^(p_25) (^4.7)-' - 1.2 X 102^(r4.7)'/'(P-25)-'A;^] « (14) 

for k ~ k{njnax) (eq-(6)) and T, p close to our scaling values, the first term dominates 
and gives Tg ~ 2.9 x 10^ yr. This initial growth time-scale is much smaller than the 
free-fall time of the condensed blobs T//,b ~ {G pb)~^^'^ and indicates that initially 
the condensations will not collapse gravitationally, but rather contract and cool 
because of the thermal instability itself. (In the presence of B the time is given by 

2 

Tg^B ^ 9 X 10^^T4.7 (p_25)~-^(l + ^) s, which for v\/Vg < 1 is also of the order of Tg.) In 
fact, the condensed blob will become gravitationally unstable if its mass exceeds the Jeans 
mass Mj = 1.2 (/c^ Tb/z^mi^)^ P~^/^ (where is the initial temperature of the blob 
and P is the pressure of the unperturbed medium), which is much larger than Muob for 
temperatures ~ 10^ K. 

A blob in pressure equilibrium with the ambient gas will have a density 
Pb — {T jib/Tbp) p (where pb ~ 1-22 for a weakly ionized gas inside the blob). 
Once the temperature of the blob dropped below 10"^ i^, the cooling will be dominated by 
the coUisional excitation of the rotational and vibrational transitions of the H2 molecule (the 
heavy elements also contribute, but for very small abundances Z/Zq they do not dominate 



-9- 



the cooling, see, e.g.. Fall & Rees 1985). Thus, just after its formation the blob will 
cool in a time rc,i ~ {?> pi,kBTij){2 jii,mH Lh^)~^ ; with Lh^ ~ 10"'^'^ nb{H2) erg cm~^ 
for an initial hydrogen number density in the blob nb{H) ~ 1 cm~^ (e.g. Lepp & 
ShuU 1983). The number density of molecular hydrogen nb{H2) can be evaluated from 
statistical equilibrium taking into account the relevant reaction rates and the effect of 
UV dissociation due to the hotter ambient gas; wc find nb{H2) ~ 3 x 10^'^ xni,{H) at 
Tb ~ 9000 i^, where x = nb{H^)/nb{H) < 1 is the fractional ionization (see, e.g.. Fall & 
Rees, eq.(18) and references therein). In this estimate we have assumed a photoionization 
rate kuv — 2.4 x 10~^^ [r/kpc)^^ appropriate for our physical conditions. This cooling 
time Tc^i ~ 2.0 x lO^a;"^ (Tfe/9000iir)7/r is now larger than the condensation formation 
time Tg (eq. 14) for x < 10"^ (e.g. Dalgarno & McCray 1972, Fall & Rees 1985), but still 
much smaller than the protogalactic free-fall time tff ~ 4 x 10^ (p/10~'^^ g cm''^)~^^'^yr, 
so that the blobs will continue to cool within the protogalactic infall gas. A cold blob 
moving relative to the ambient gas may be subject to the Kelvin-Helmholtz (K-H) 
instability at the interface between blob and ambient gas with a characteristic growth time 
tk-h — \{Pb/ pY^'^ /u, where is the radius of the blob and u is the relative velocity 
(Murray et al. 1993). Buoyancy effects may accelerate some blobs to speeds comparable 
to the sound speed, for which tk-h ~ ^b/ {l ksTb/ jibf^HY^"^ ■ Imposing the "survival" 
condition for a blob against K-H instability Tc,i < tk-h gives a lower limit to the size of 
the blobs Af, > 6.3 x lO^^x-^ {Tb/'^mm f^ cm. Thus, blobs moving with relative speed 
close to the sound speed are expected to break up since they do not satisfy this condition 
for X < 10""^ (eqs. 6, 10, 11). On the other hand, those "moving with" the collapsing 
protocloud will be stable to the K-H instability. The presence of magnetic fields will further 
inhibit the K-H instability effects, although the surviving fraction of blobs must not be 
expected to be large since in that case the observed MACHO fraction would also be large 
which is not the case (Alcock et al. 1995, Gates et al.l995). A detailed calculation of the 
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K-H effects is out of the scope of this work and will be discussed elsewhere. 

Finally, the blobs must survive to a much "dangerous" later epoch in 
which the protocloud has been heated to T ~ lO^K (for some radius) and 
p ~ 10~'^'^ g cm~^ (e.g. Fall & Rees 1985). The blobs which should have by now 
cooled to Th < 1000 since Tcj — 13 nb{H) / nb{H2) yr (where we have assumed 

^ 5 X 10~^'^ nb{H2) erg cm.~^ for nb{H) ~ 10^ cm"^; see Lepp & ShuU 
1983) will be heated by X-rays from the hot ambient protocloud. The effective 
cross section for photoionization of the blob gas due to ~ 0.1 KeV photons is 
a ~ 1.8 X 10~^° {hu / 150 eV)~^ cm^ . Since the radiation is mainly thermal bremsstrahlung 
with an intensity oc exp{—hv/kBT), about 90% of the heat will be deposited in a layer of 
width /x ^ (nfc(/J)(T)-i ~ 1.3 x 10^^ {T / 10^ Kf{p/ 10""^^ g cm^^y^Tb/ 1000 K) cm 
by photons having < /iz/ > ~ 2.3A;bT. This length scale defines a mass 
Mx ^ 2.1 X 10-^ {T/ 10^ Ky{p/ 10-"^ g cm-^)-'^ {nl lOOOKf Mq. Therefore, blobs 
with Muoh < Mx will be heated throughout and the more massive ones only in a thin 
layer ~ Ix near their surfaces. This X-ray heating will kill the smallest blobs from the 
original spectrum (eqs. 8 and 12), providing a possible explanation to the lack of low- mass 
(M < 10-^ Mq) candidate events for MACHOs (Zylberajch 1994). These figures must be 
considered as lower hmits to the heating rates since any other sources such as supernova 
explosions or cosmic rays would increase Ix- However, we must stress that the blobs are in 
principle able to cool quickly to < lOOK in a time-scale ~ 3.3 x 10^ [nb{H) / nb{H2)Y yr (for 
Lh^ ~ 2 X 10"^^ nb{H2)'^ erg cm''^ for nh{H) ~ 10^ cm"^; see Capuzzo-Dolcetta et al. 
1995), which is shorter than the protogalactic t// at that epoch. These low temperatures 
are enough for the blobs to become gravitationally unstable, and therefore it is not unlikely 
that a number of them can survive until today. 

Several other processes (not discussed here) may be also important for the evolution of 
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the blobs (including ambipolar diffusion at < 10^ K (e.g. Shu et al. 1987), and ~ 10^ 
collisions between blobs since their formation), which may bias the final mass distribution 
in a yet unknown manner. 

The main result of these calculations is that a thermally unstable regime of the 
protogalactic gas before the standard stellar formation trigger induces the formation 
of condensations or blobs in a suitable mass interval. If the physical conditions of the 
protogalactic gas are not very different from the standard models of the protogalaxy, 
the most abundant surviving blobs have typical masses M > IQ-^ Mq, consistent with 
the microlensing candidates reported by the MACHO and EROS collaborations. In this 
scenario, the MACHOs may be identified with cold clouds which started to collapse from 
condensed blobs because of the thermal instability of the gas. It has also imphcations 
for the generation of low-mass Pop III luminous stars, a subject that has received some 
attention in the recent past (e.g. Silk 1983, Bessel & Norris 1984 ; Beers et al. 1985 ; Gass 
et al. 1988). 

After the completion of this work we had notice of a recent work (De Paolis et al. 
1994) which discusses the issue of thermal instabihties forming M ~ 10^ hot clouds 
later fragmenting into low-mass objects identified with the MACHOs. Two comments are 
in order. First, as stated by the authors, clustering of the lensing objects is expected and 
this prediction can be checked against accumulated observed events. Second, the clustered 
MACHOs formed inside those massive clouds should in any case add up to the fraction 
resulting from the previous epoch as discussed here; thus we may be observing different 
dark populations originated closely in time but by different physics. 

We are indebted with the referee A. Ferrara for his fruitful comments and suggestions. 
Discussions with R.Opher, S. Viegas, A.C.S. Friaga, A. Dal Pino Jr., and G.A. Medina- 
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